Fitness in time-dependent environments includes a geometric contribution 
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Phenotypic evolution implies sequential fixations of new genomic sequences. The speed at which 
these mutations fixate depends, in part, on the relative fitness (selection coefficient) of the mutant 
vs. the ancestor. Using a simple population dynamics model we show that the relative fitness in 
dynamical environments is not equal to the fitness averaged over individual environments. Instead it 
includes a term that explicitly depends on the sequence of the environments. This term is geometric 
in nature and depends only on the oriented area enclosed by the trajectory taken by the system in the 
environment state space. It is related to the well-studied geometric phases in classical and quantum 
physical systems. We discuss possible biological implications of these observations, focusing on 
evolution of novel metabolic or stress-resistant functions. 
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I. INTRODUCTION 

Organisms react to long-term changes in environmen- 
tal conditions by sequential fixation of new genome 
sequences, mostly corresponding to increasingly more 
adapted phenotypes. However, often environmental 
changes are faster than the characteristic time for 
mutation-selection cycles needed to evolve an optimal 
phenotype. In such cases, depending on the structure and 
time scales of the fluctuations, a dynamic environment 
creates dynamic fitness landscapes [T], promotes sensing 
[2] , modularity H] , switching [5] , and can change the 
speed of adaptation [BJ [7] . 

The effect of fluctuating selection and/or population 
size on the population-genetics dynamics have been ex- 
tensively studied over the years [7j [5] , starting with the 
introduction of the concept of adaptive topography by 
Wright [|k More recently, the evolutionary dynamics 
of density regulated populations in fluctuating environ- 
ments has been elucidated in ecologically realistic mod- 
els [TDHH]. These bridge the gap between the classical 
population dynamics exhibiting very diverse responses to 
fluctuating environments [T3l [14] and classical population 
genetics models. However, a complete understanding of 
the effect of fluctuations on population and evolutionary 
dynamics has not been achieved yet. 

Some of the relevant parameters describing evolution- 
ary response of a population to a changing environment 
are the rate at which new genotypes are created (muta- 
tion rate), the relative fitness of new phenotypes, and the 
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total population size. We concentrate on the case of envi- 
ronments changing on scales longer than an individual's 
lifetimes. This is relevant, in particular, for bacterial 
populations confronted with daily environmental changes 
(natural or artificial) [15] . for longer- living organism af- 
fected by seasonal variations, or for pathogens experienc- 
ing transmission, uncontrolled growth in a new host, and 
then effects of the host immune system. For example, 
in the now-classic long-term E. coli evolution experiment 
[IB] , bacterial cultures are diluted daily, and the environ- 
ment (i.e., cell growth and death rates) changes during 
dilution events and between them due to depletion of re- 
sources, cell density growth, and cell-to-cell interactions. 
These experiments are a great model to study clonal com- 
petition [17]. Interestingly, the number of accumulated 
beneficial mutations is relatively small, considering that 
every single point and many possible double mutations 
have happened thousands of times in the 25-year history 
of the experiment. This discrepancy is likely largely ac- 
counted for by strong bottlenecks at dilution times, when 
most new mutations disappear by chance. However, all 
clones, even beneficial ones, experience additional huge 
fluctuations in their reproductive rates during the course 
of the experiment. It remains to be seen if such fluctu- 
ations can contribute to the slowing down of the evolu- 
tionary adaptation as well. 

In this article, we make a step in this direction by 
studying effects of fluctuating environments (represented 
by birth and death rates) on the effective selection coef- 
ficient. Using analytic and computational tools, we in- 
vestigate a model of a heterogeneous population (a back- 
ground strain and a newly emergent mutant) under the 
assumption that the time scales of the clonal frequency 
dynamics on the one hand and the environment fluctua- 
tion on the other are both much larger than the division 



2 



time, but not necessary well separated from each other. 
We start by showing that the selection coefficient in an in- 
finitely slowly changing environments is given by a time- 
average of static selection coefficients corresponding to 
each environment. However, for environments varying 
at a slow but finite rate, such timc-avcragc is not the 
whole story. A new contribution emerges. For exam- 
ple, in a cyclically oscillating environment, this contri- 
bution to the selection coefficient is independent of the 
speed of variation and depends only on the sequence of 
environments visited during each cycle. The contribu- 
tion is non-zero only for nontrivial coupling between the 
environment and the population dynamics, represented 
as a multi-dimensional trajectory in the space of birth 
and death rates. The contribution changes sign when 
the sequence of the visited environments is reversed, ft 
is largely independent on the speed of the dynamics. Fi- 
nally, it scales quadratically with the amplitude of the 
environmental fluctuations. In other words, the contri- 
bution is geometric in nature. We believe that this has 
not been noticed before in the context of population dy- 
namics. 

We will focus on the deterministic approximation to 
population dynamics. Geometric effects are well-known 
for slowly changing deterministic dynamical systems [181 
I19j . While evolutionary dynamics of a population driven 
together by forces of mutation, drift and selection cannot 
be accurately described deterministically, we believe that 
our model is meaningful even for a stochastic case for 
large population, low mutation rate, and strong selection. 
Indeed, the recent observation that stochastic dynamical 
systems are also subject to geometric corrections suggests 
that deterministic vs. stochastic treatment of population 
dynamics is not crucial for the phenomenon [20l |2T] . 

In what follows, we develop our results in a relatively 
simple two species population model with bilinear, sym- 
metric competition, which we believe is general enough 
to capture the main effects of fluctuations for a large class 
of related models. We first solve the system in the limit 
of small differences between the birth and the death rates 
of the competing species. We derive expressions for the 
selection coefficient in the limit of stationary, very slowly 
continuously, and infrequently discontinuously varying 
environments. The selection coefficient for arbitrary time 
scale of the environment fluctuations can be derived then 
using a perturbative approach. 



II. MODEL 

Let Xi be the number of individuals of genotypes 
i, i = 1,2, in a large asexual population. We assume 
that Xi ^> 1, so that demographic (phenotypic) fluctu- 
ations and random genetic drift can be neglected. Wc 
refer to x\ as an ancestral phenotype, and to x 2 as a mu- 
tant. The competition between the two is described by 
a driven two-dimensional Lotka-Volterra (logistic) model 



[HE] 



±i = x\ [&i(t) - d\[t){x\ + x 2 )] , 
x 2 = x 2 [b 2 (t) - d 2 (t)(xi +x 2 )}. 



(1) 



Here bi(t) represents the birth rates, and di(t) parameter- 
ize the death rates for each of the genotypes. Generally, 
all parameters are time dependent. 

Following classical models of ecological population ge- 
netics, we view our model as a particular form of the 
more general dynamics. Defining the total population 
size, x(t) — x\(t) + x 2 (t), we write 



i'2 



xigi(x,t), 
x 2 g 2 (x,t). 



(2) 



Here g is the generalized growth rate. For this system 
of equations to represent the dynamics of a realistic self- 
sustaining population, gi{x) must be negative for large x, 
and it must have at least one zero. Our approach applies 
to a very general subset of such growth rate functions 
provided that the system, Eq. pi), has exactly one fixed 
point on each of the axes x% = in addition to the trivial 
unstable extinction point (0,0). 
One traditionally takes [H] 



9i 



i(*) 



f 



/(*) 



f(Ki(t)) 



(3) 



where 's are the intrinsic maximum growth rates of each 
genotype, if unconstrained by limited resources. The 
terms r i f(x)/f(K i (t)) represent the reduction of these 
rates due to competition for resources. This reduction 
depends only on the total population size x(t) and on K%, 
which are stable total populations of the isolated pheno- 
types i supported by stationary resource-limited environ- 
ments. K' s are referred to as the carrying capacities. Our 
approach applies for any non-negative, monotonously in- 
creasing f(x), as explained above. However, for sim- 
plicity, we now concentrate on f(x) — x. In this case, 
the competition is linear and symmetric, and the simple 
Lotka-Volterra model is recovered with 6j(t) = r,(t) 
and d l {t)=r l {t)/K l {t). 

We are interested in modeling competition of the an- 
cestral genotype with the mutant one. The two are very 
close in the genotype space, essentially one mutation 
away. Since mutation effects are, in general, small |17j . 
we assume that the differences between g\ and g 2 are also 
small, 



giQM) -ff2pM) 
9i(x,t) + g 2 (x,t) 



< € < 1. 



(4) 



This corresponds to small differences in the parameters 
bi, di, Ti, Ki. We assume this from now on. In particular, 
it is possible that differences between the mutant and 
the ancestor parameters at any particular time are much 
smaller than the variations of each of the parameters over 
time. 
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III. PRELIMINARIES 

In order to determine the conditions under which the 
mutant, initially present in small numbers relative to the 
ancestor, invades the population, we explicitly integrate 
the model, Eq. Q. We write the dynamics of the total 
population size x = x\ + x 2 : 



b\(t)xx b 2 (t)x 2 



di(t)xi d 2 (t)x 2 



(5) 

To the zeroth order in e < 1, this does not depend on 
the individual values x\ and x 2 . 



x = x [b(t) - d(t)x] + 0(e), 

where we have defined 

bx(t) + b 2 {t) 



b(t) 



We also define 



d(t) 



d l {t)+d 2 {t) 



P 



x 2 



Xi+X 2 



(6) 



(7) 



(8) 



the fraction of the mutant in the whole population. This 
obeys 

P = p(l- P ) {[6a (i) - hit)} - [d 2 (t) - d 1 (t)] x} . 
The model then reduces to 



x = x [b(t) — d{t)x] , 

p = p(l-p)[5b(t)-Sd(t)x] 



(9) 



where we have used the notation 5(b, d) for small (order 
e) time dependent differences between the corresponding 
mutant and ancestral rates. To simplify the notation, for 
any pair of parameters (Vi,V 2 ) describing the ancestor 
and the mutant, we write V = (Vi + V 2 ) /2, and 6V = 
V 2 —V\. In addition we always assume \SV/V\ = O(e) <C 
1. 

To the zeroth order in e, the dynamics of the total 
population size defined by Eqs. (JlJ is now uncoupled from 
the dynamics of the mutant fraction 



x(t) 



x(0)e^ dTb ^ 
x(0) J^dt'd(t')efa' dTb ^ 



(10) 



Due to the small variation assumption, Eq. Q, p(t) 
changes on time scales much longer than x(t). On these 
time scales, x(t) converges to a unique (up to the first 
order in e) attractor x a (t), independent of the initial con- 
ditions, 



x a {t) 



(11) 



J^JLt' d(t<) e£ dTb(T) 
Then the slower dynamics of p is 



logitp(t) = logitp(O) + / dr [Sb(r) - M(t):e(t)] , (12) 
Jo 



where logitp = \ogp — log(l — p). The obvious first les- 
son from this equation is that the clone with the largest 
average growth rate, (<?j) > ^ dt [bi(t) — di(t)x(t)} for 
some large T, will have an advantage. 



IV. SELECTION COEFFICIENT 

For coefficients varying periodically with a period T, 
we write for the logarithmic change of the mutant-to- 
ancestor ratio, logitp, over time T ^> T, 



A(T) = logitj>(T)-logitp(0) 



T 



fidr [Sb(r) - 5d(r)x(r)} 



T 



sT, (13) 



where the last equality defines the selection coefficient, 
s. It is the sign of s that decides the stability of the 
fixed points p — 1 and p — 0. For example, for s > 0, 
p = is unstable, and the mutant phenotype invades the 
population towards a stable fixed point p = 1. 

In a constant environment, and for e -C 1, the selection 
coefficient s can be rewritten in terms of the ecological 
parameters defined in Eq. ([3| 



SK 
~K' 



(14) 



We have a classical result that selection favors pheno- 
types with larger carrying capacities (larger Ki) inde- 
pendent of the magnitude of the intrinsic growth rates rj 
PHI E]. To derive this, we rely on the fact that the total 
population is given at all times by K, and it is indepen- 
dent of the frequency of the mutants in the population. 

In this paper, we are interested in the values of the se- 
lection coefficient for temporally varying environments. 
As a consequence, the selection coefficient is now given 
by the interaction between several varying quantities. To 
simplify the discussion, we focus on limiting cases of large 
time scale separation between the environment fluctua- 
tions and individual lifetimes. 

In the regime of infinitely fast environmental fluctua- 
tions, for T — > oo, we approximate the general driven 
model, Eq. as 



x 2 



xi{gi{x)) T , 
x 2 (g 2 (x)) T - 



(15) 



We assume here that the environment variation attains 
a well defined, constant average for every state (xx,x 2 ). 
We denote this by (■ ■ -)t, where the subscript T stands 
for averaging over a period. We assume that x does not 
change appreciably over this time. For the specific case 
of the Lotka-Volterra model, the selection coefficient for 
fast fluctuations, Sf, can be computed using the formula 
for the constant case, Eq. (14), keeping in mind that one 



has to use the average values of the relevant coefficients: 



s f = 6(r) -S(r/K) 



(r) 



\r/KY 



(16) 
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In the opposite limit of an infinitely slow parameter 
variation, the total population is equal to the carrying 
capacity at all times, x(t) = K(t). In this case, the 
quasi-stationary (qst) selection coefficient s is 



fqst 



T , , JK(t) 
dtr ^KW 



dts(t), 



(17) 



where the period T is much longer than the individual's 
lifetime. This allows for a proper average to be attained. 

In both limits, the sign of the selection coefficient does 
not depend on the average carrying capacity [H] H3 • In- 
deed, it is possible to have a slowly varying environment, 
in which the mutant has, on average, a larger carrying 
capacity but a lower fitness. In both limits, the selection 
coefficient becomes independent of the speed of environ- 
mental variations, and it is symmetric with respect to 
time reversal for the driving parameters. 



V. CONTINUOUS, DETERMINISTIC, 
OSCILLATORY ENVIRONMENTS 

We now proceed to a more realistic case of an environ- 
ment fluctuating slowly, but not infinitely slowly, com- 
pared to an individual's lifetime. This condition allows 
us to derive a perturbative approximation for the selec- 
tion coefficient valid when b(t),r(t) ^> A are satisfied at 
every t. Our approximation is based on a simplified solu- 
tion for the dynamics of the total population size x a (t), 
Eq. (If). By making a variable change y(t) = J*dTb(r), 



we write 



1 



1 



rV 
J — c 



3 b(z) 



I. 



oo K(z 



(18) 



In the limit of slow environmental changes, the carrying 
capacity K(y) varies slowly, and the integral in the de- 
nominator is dominated by the value of l/K(z) around 
z = y. In this regime, 



f 



1 



KM 

K(z) K(y) K*(y) 



(z - y) for (y-z)< y. (19) 



Using Eq. (19 1, we now derive an approximation for 



the total population trajectory x a valid in the qst regime. 
We denote it as x oa , 



x qa (£) ~ K(t) 



K'it) 
r(i) ' 



(20) 



This solution represents the correction to the quasi- 
stationary result x qs t(£) = K(t) as a first order perturba- 
tion in the small ratio between the rate of change of the 
environment and the typical rate of change of the total 
population. Note that the approximation is consistent 
with the intuition that the instantaneous total popula- 
tion falls behind the instantaneous carrying capacity. 



The selection coefficient can be expressed now as 



'qst 



where 



jgeom 



1 

T 



dr 



5r(r) 5K(t) 



r(r) K{t) 



K'{t) 
K(t) 



(21) 



(22) 



is a geometric contribution to the selection rate. The 
geometric nature of this term can be better understood 
if we express the change in the mutant-to-ancestor ratio 
as 



A(T) 



'qst 



T+A E 



>(T). 



(23) 



We note that, for any reparameterization of time, A = 
X(t), Ageom can be written in a very similar form 



Ag CO m (7") 



A(T) 



5rjX) _ 5K{\) 

r(x) k(x) 



K'(X) 
K(X)- 



(24) 



which emphasizes that it depends on the trajectory it- 
self, ri 2(A), Ki,2(X), rather than on how this trajec- 
tory is traversed. As any closed contour integral expres- 
sion, this expression can be transformed into a surface 
integral over any 2D domain bounded by the trajectory 
[ri(i), r2(t), Ki(t), K2(t)] in the parameter space. In par- 
ticular, using variables 



X = Slog 



K 



y = log k 



(25) 



and the Stokes theorem, we can equate A geom (T) with 
the oriented area bounded by the trajectory for times 
t G (0,T) in the plane (X,y). 

In other words, A goom is a truly geometric term in the 
spirit of geometric phases in quantum or classical me- 
chanics [HI [19] . The geometric nature of the change in 
the population composition over long times, Eq. (24), is 
the main result of the paper. It allows us to make impor- 
tant macroscopic predictions about the population dy- 
namics that will hold generally irrespective of the micro- 
scopic details of the model. First, the geometric changes 
in the relative fraction of the mutant depend on the se- 
quence of the environmental states in addition to their 
identity: same environmental states may have very dif- 
ferent effects depending on the order in which the states 
are visited. At an extreme, a reversal of the order (time- 
reversal) would change the sign of the geometric contri- 
bution, which may make a deleterious mutation advan- 
tageous, and vice versa. To our knowledge, such depen- 
dence of the effective selection coefficient on the sequence 
of the environmental states has not been noticed before 
in population biology. Second, the contribution to A geom 
depends only on the oriented area covered in the param- 
eter space (and thus, in particular, on the number of 
periodic oscillations), but not on the speed of traversal 
of the trajectory. Figure [T] illustrates these features: even 
when the environmental dynamics involves backtracking, 
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FIG. 1: Mutant fraction as a function of time for two 
sample environment trajectories (X(t),y(t)): X\(t) = 
0.02 sin [uit + sin(wit)] , yi(t) = 1 + 0.1 cos [uit + sin(wit)] 
and X 2 (t) = 0.02 sin [u 2 t + 2.5 sin(cj 2 t)] , ya[t) = 1 + 
O.lsin^t + 2.5sin(o;2t)] where wi/1.4 = o; 2 = 27r/40. (A) 
The two trajectories for 3^ = log K are shown; the first has 
the frequency 1.4 times the second, and the second reverses 
twice before completing the full cycle. (B) Nonetheless, the 
shapes of the trajectories (X(t),y(t)) are the same for both 
examples. (C) Instantaneous and one-period-averaged mu- 
tant fractions for both trajectories. The average growth of A, 
given completely by a geometric term, is linear. The slopes of 
the two curves are different by exactly 1.4, so that A is only 
dependent on the number of elapsed periods. This is indicated 
by the horizontal line connecting the two averages delayed by 
the same number of periods. Thus the geometric contribu- 
tion to the mutant fraction depends only on the shape of the 
contour in the parameter space and on the number of cycles, 
but is independent of the speed of the trajectory traversal. 
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FIG. 2: Illustration of the geometric nature of the mu- 
tant fraction dynamics. (A) Three different trajectories 
(X(t),y(t)): Xi{t) = -0.02cos(wt), yi{t) = l + 0.1sin(urf), 
X 2 (t) = O.Olcos(wt), y 2 (t) = 1 + .lsin(wt), X s (t) = 
0.02cos(W), y 3 (t) = 1 + 0.1sin(2wt) where u = 2tt/100. 
The second trajectory (solid line) encloses exactly half the 
area of the first (dotted line), and the two are traversed in 
opposite directions. The oriented area enclosed by the third 
trajectory (dash-dotted) is zero. (B) The average mutant frac- 
tion change for the first trajectory is equal to the oriented area 
and is, therefore, twice that for the second one, and in the op- 
posing direction. The quantity is zero for the third trajectory. 



The state occupied between t a and t a +i will be denoted 
by \i a . We assume that the interval (f a +i — ta) is l° n g 
enough so that the total population x{t) reaches the 
carrying capacity long before the environment switches 
again, that is l/r Mi <C (t a+1 — t a ). In this case one can 
derive the qst contribution as a sum over all of the envi- 
ronment states 



^qst 



t a ). 



(26) 



At each switch, there is an extra contribution because 
x(t > t a ) reaches the value K^ a with a delay. That is, 



from Eq. (10), we derive: 



the overall contribution per period still does not change. 
The dependence on the area in the parameter space also 
suggests that the geometric contribution scales as the 
square of the fluctuation amplitudes. Finally, to achieve 
a nonzero area, more than one parameter must be chang- 
ing, and they must change incoherently. We illustrate 
some of these features in Fig. [2] 



VI. SWITCHING AMONG DISCRETE 
ENVIRONMENT STATES 



The approach can be extended to a more com- 
mon model of piecewise constant environments, see 
e.g., Refs. [21 U\- Consider the case of parameters 
abruptly changing between m sets indexed by fx = 
1 . . . to, (r^, r% , K± , K% ), at possibly random times t a . 



x(t a < t < t a+1 ) = 



KHa-1 



3 r"-(f-t) +1 



(27) 

Integrating Eq. (27 1 results in a geometric contribution 



after M environment state changes 



A 



gOOlll 



Sr^(X) 
r^(A) 



SK^(X) 



log 



K* a (A) " 

^a-l(A) 

(28) 

The fact that Eq. ( 28 ) is independent of the actual time 
spent in each state and depends only on the sequence 
of environmental states is the signature of its geometric 
nature, illustrated in Fig. [3j Importantly, unlike in the 
continuos variation case, Eq. p4[ ), A goom in Eq. (28) can 
have a finite value even if parameters change only be- 
tween two states. Hence it is unclear if the contribution 
can be interpreted as an oriented area enclosed by the 
trajectory in the parameter space. 
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FIG. 3: Mutant-ancestor competition for environment fluctu- 
ating between two states. (A) Time dependence of x\ and x 2 - 
(B) Time dependence of the logarithm of the population ratio, 
A(T). The two states are characterized by (n = 2.3, r 2 = 
2, K X =K 2 = 1) and (n = 1.1, r 2 = I, Kx = K 2 = 2). The 
time spent in each state is uniformly distributed between 4 
and 10. The equal carrying capacities ensure zero qst contri- 
bution while the mutant x 2 is winning in the long-term due 
to the geometric contribution. The evolutionary pressure is 
exerted only at the very beginning of the residence period 
in each environmental state, and the total population and 
the mutant fraction stay constant for the rest of each phase. 
Thus the mutant ratio drift depends only on the number of 
switches, but not on the duration of the process. 



VII. CONTINUOUS STOCHASTIC 
ENVIRONMENTS 



Now using suitable continuity properties of the parame- 
ters' trajectory, we transform the geometric contribution 
to 



i gcom 



(T) = 



T 



dt 



A A 

E E 

a=l f)=a+l 



(31) 

are clear from 
Eq. ([311: A gcom (T) depends only on the length of the 



The geometric properties of A g 



parameters' trajectory, is antisymmetric with respect to 
time reversals, and is nonzero only if multiple param- 
eter vary simultaneously and incoherently. Note that 
Eq. ( |3ip is valid only for parameter variations with small 
(bounded) speeds. Therefore, if the parameter dynam- 
ics, 7 a (i), are modeled as multidimensional Wiener pro- 
cesses, care must be taken to regularize and properly de- 
fine the stochastic integrals in Eqs. (29 31) j2"2"] . 

Equations (29) and (31) represent a natural extension 



of the geometric correction to acyclic trajectories (23j . 
While now the geometric term A gcom (T) is aperiodic, for 
parameters dynamics with a stationary distribution of j a 
and 7 a , A ge om(T) still has a mean linear dependence on 
T for large times: 



lim Agcom 



(T) 



r 



A A 

E E > 

a=l /3=a+l 



dC aP (t) 



dt 



(32) 

where C a p{t) = (7«(0)7 l g(t)) are time dependent corre- 
lations of the environment. Note that the derivatives 
dc a ji{t) inversely proportional to the correlation 

t=o 

times of the process. Moreover one can identify the terms 
in the rhs of Eq. (32 1 as products of the Berry curvature, 
Kp a — K a p, previously introduced in the classical and 
quantum geometric phases literature [HI HI], and, for 
a f3, the rates of growth of the oriented areas bounded 

t=0 



by the process 



dt 



Often environments change in a continuous but unpre- 
dictable way, such that the typical rate of change is still 
small. This scenario is modeled with Gaussian fluctua- 
tions of the parameters [HI [52] . Denoting all parameters 
with a single symbol j a , a = 1, . . . , A, we generalize our 
result, Eq. (24), and represent the geometric contribution 
for randomly driven Eqs. |2]) as a line integral [2TJI |2"TI |2"3"] 



x(T) 



r 



dt 



A 

E 

a=l 



fa(jl(t),...^ A (t))j a (t). (29) 



Here T is a long time that allows for averaging, and f a are 
some model-dependent functions. Since fluctuations are 
small, we expand f a to the first order in the fluctuations 
around the average parameters 

A 

K a f3jp(t)- (30) 

0=1 



VIII. POSSIBLE EXPERIMENTAL EFFECTS 

The existence of geometric corrections to fitness in a 
time dependent environment requires that changes in the 
environment are felt by the population on multiple time 
scales. In the model, Eq. Q, the immediate change in 
the growth rates and the delayed effect of the population 
reaching the carrying capacity provide these scales, but 
other mechanism would work as well. Similar effects will 
be encountered in almost any situation when a popula- 
tion responds to asynchronous changes in multiple exter- 
nal stresses or nutrient supplies. Therefore, the geometric 
effects must be considered when modeling emergence or 
fixation of new metabolic or stress-resistance functions in 
the presence of environmental changes. We suggest that 
the relative timing of fluctuations of extracellular nutri- 
ent /stressor concentrations will affect the relative fitness 
advantage of these functions. 
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FIG. 4: Simulated dynamics of the logarithm of the relative 
population size for two partially antibiotic resistant popula- 
tions competing for the same consumable carbon source in 
a chemostat. The concentration of the antibiotic changes 
as A(t) = 0.1 + O.lcos(ait) in arbitrary units. The nu- 
trient influx is 1.0 + cos(ujt + <j>). The nutrient is cleared 
by the chemostat and consumed by both strains in propor- 
tion to the population growth, resulting in the concentration 
v(t). The growth rate of either population is proportional to 
V (1 + A/Ka)~ 1 (1 + Kv/v)~ l . V is the maximum growth 
rate, K v is the Monod growth constant, and Ka is related to 
the minimal inhibitory concentration for the antibiotic. Ka 
for the more resistant strain is 14% higher than for the less 
resistant one, but its V is 5% smaller to account for the cost 
of resistance |24) . The numbers are chosen such that the av- 
erage growth for very slow environmental changes (solid line) 
is almost the same for both strains. Depending on the phase 
4> (dashed and dash-dotted lines), either the resistant or the 
non-resistant strain has the higher growth rate and will even- 
tually take over the population. The "difference" line shows 
the nearly linear difference between the strain fractions for 
the two opposite environmental trajectories. 



Of a particular interest is emergence of antibiotic resis- 
tance in bacteria. Mutations conferring antibiotic resis- 
tance often decrease ability of cells to grow in the absence 
of antibiotics, but provide a growth advantage in their 
presence [21] • At the same time, delivery of antibiotics 
is hardly ever uniform, and nutrient supplies also fluc- 
tuate. Focusing for simplicity on periodic nutrient and 
antibiotics concentration changes, we see that the time 
delay, or the phase lag, between the changing concentra- 
tions will join their amplitudes and the period in selecting 
whether a resistant strain will fix or not. We illustrate 
this in Fig.[4j depending on the phase difference between 
the nutrient and the antibiotic influx, either the resistant 
or the faster growing bacterium will be selected for. A 
robust prediction of our theory is that the difference in 
the logarithmic fractional population changes between an 
environmental trajectory and its time reversed counter- 
part will grow almost linearly in time with the number 
of periods. We emphasize that the effect is different from 
episodic selection [25] , where only frequencies and mag- 
nitudes of antibiotic selection episodes determine fixation 
of the resistant strain. 

Another experimental system where our predictions 



can be important is evolution of a metabolic pathway 
corresponding to a new metabolite, when both the old 
and the new metabolite concentrations change in time. 
In such a case, one would need to take into account pos- 
sible effects of catabolite repression and di-auxic growth 
in addition to instantaneous effects of metabolite concen- 
trations on the birth/death rates. Nevertheless we expect 
that careful modeling of these effects will also uncover the 
fitness sensitivity to the timing of pathway activation. 

An important characteristic of the geometric effect is 
that it is much harder to be observed in typical serial 
dilution experiments, especially when the environment 
changes are only imposed at the dilution points. Such 
experimental protocols will miss important effects that 
may be relevant for wild-type conditions. 



IX. DISCUSSIONS 

Fixation dynamics of mutants in a large class of math- 
ematical models is governed by a single effective param- 
eter, the selection rate, obtained as a time average of the 
instantaneous growth rate difference between the mutant 
and the ancestral population. In population dynamics 
with symmetric competition, and in the limit of small dif- 
ferences between the mutant and the ancestor, the total 
population size is decoupled from changes in the popula- 
tion composition. Instead the total population enters the 
fixation dynamics only as a time dependent parameter. 
Then the population growth rates and the selection coef- 
ficient depend on the interplay between the time scales of 
the population dynamics and the environmental fluctua- 
tions. For infinite separation between the time scales, the 
selection depends only on values of environmental param- 
eters. More specifically here the fitness difference can be 
expressed as a function of growth rates and carrying ca- 
pacities averaged over all of the environmental states and 
independent of the period of the fluctuations. Nonethe- 
less, due to the non-linear dependence of the growth rates 
on the environmental parameters, the average fitness dif- 
ference is not necessarily the same as the fitness difference 
for the average environment. 

This quasi-steady state approximation breaks down for 
faster environmental changes. The mutant fraction dy- 
namics is now dependent not only of the period of en- 
vironmental changes but also on the sequence of succes- 
sive environmental states. In particular the first non- 
adiabatic correction is always anti-symmetric with re- 
spect to time reversals, and it is geometric in nature. 
As long as the fluctuations in the parameters are large, 
this non-adiabatic correction can be of the same order 
of magnitude in the birth and death rates variation as 
the qst contribution to the fitness difference. The geo- 
metric nature of this term constrains the effect that en- 
vironment fluctuations can have on fitness differences. 
Indeed, as other geometric contributions [18j [19], this 
effect is independent of the instantaneous speed of varia- 
tion of parameters. In ecological terms, this implies that 



the geometric contribution to the mutant ratio drift does 
not depend on how fast the environment changes, but 
only on the sequence of environmental states. We illus- 
trate this in Figs.[TJ[3j Further, we note that the mutant 
fraction drift, A, can be seen as a line integral in the 
parameter space, cf. Eq. (24). This implies that only 



multidimensional and off-phase parameter variations can 
give nonzero long-term contributions to the population 
dynamics. 

For the results derived in this work, the assumption 
of an oscillatory environment is not essential. Our con- 
clusions, and the concept of geometric phase in general, 
are valid for non-cyclic environment dynamics |23j . Typ- 
ically such dynamics is represented with a Gaussian and, 
in general, uncorrelated noise [51[TTJ[T2]. While a detailed 
extension of the present results to random trajectories is 
beyond the scope of this paper, we have shown here that 
the geometric contribution to the selection coefficient is 
present generically if and only if the population dynam- 
ics contains multiple correlated parameters driven by a 
colored noise, cf. Eq. (32). 



In this article we have focused on deterministic popu- 
lation dynamics with small parameter differences among 
the competing species, which is equivalent to frequency 
independent selection. We expect that a similar geomet- 
ric phase contribution to the fixation dynamics is present 
in stochastic Fisher- Wright type models, as well as mod- 
els that exhibit various frequency and density dependent 
selection effects. 



The results in Eq. (24) allow us to make a conclusion 



that is independent of the exact variation of the param- 



eters and the exact details of the model. Namely, for 
clones with the same mean fitness, the clone that has 
a higher growth rate when the environment is abundant 
(increasing carrying capacity) will have a selective ad- 
vantage over the clone that performs well when the car- 
rying capacity decreases. This is important during acqui- 
sition of new metabolic or stress-response functions, as 
discussed above. Further, in the case of the long-term E. 
coli evolution experiment [T5], we point out that unless 
mutations manifest themselves in a positive way during 
the exponential growth phase following a serial dilution, 
daily variability of the environment would make it harder 
for mutations to fixate even without stochastic effects as- 
sociated with the dilution bottlenecks. 

We conclude with an observation that species with the 
fitness advantage in the average environment, with the 
average fitness advantage over all environments, and with 
the average fitness advantage for a particular time course 
of the environment are not necessarily the same species. 
In particular, a naively deleterious mutation can fixate in 
a population due to these temporal effects. We believe 
this to hold true independently of many of the simplifying 
assumptions of our toy model. 
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